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A reformulation of the Tliirring model as a gauge theory on both continuum space- 
time and discretized lattice is reviewed. In (1+1) dimensions, our result reproduces 
consistently the bosonization of the massless Thirring model. In (2+1) dimensions, 
the analysis by use of Schwinger-Dyson equation is shown to exhibit dynamical 
fermion mass generation when the number N of four-component fermions is less 
than the critical value A^cr = 128/37r^. 



1 Gauge Theory Formulation of Thirring Model 

The Lagrange density of Thirring model is given by 

a a,b 

Here ipa is treated as a formal one in arbitrary dimensions, but a four-component 
Dirac spinor in (2-1-1) dimensions and a, b are summed over from 1 to N. (For 
the notations and the detailed calculations, see Ref. [1,2].) Since the form of 
four- fermion interaction term is a contact term between vector currents,^ a 
well-known technique to facilitate 1 /A^-expansion is to rewrite Eq. (|l|) by intro- 
ducing an auxiliary vector field A'^ such as 

>Caux = V ^all^id^ - + ^A^AI^. (2) 

Gauge-noninvariant as the above Lagrangian is, however one may be tempted 
to regard as a gauge field. A systematic way to construct the U{1) gauge 
theory, but is gauge equivalent to Eq.(0) is to elicit the fictitious Goldstone 
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degree (or cquiYajkjoily the Stiickelberg field) based on the principle of hidden 
local symmetrylj'QEllj: 

-ChLS = ^at^D^^a + ^(^P - ^d^^^f- (3) 

a 

It is obvious that Eq.(^) possesses a U{1) gauge symmetry and the gauge- 
fixed (unitary gauge) form of it exactly coincides with Eq.(j2|), so does the 
original Thirring model in Eq.(^. The gauge theory formulation of Thirring 
model in Eq.(^) looks like quantum electrodynamics (QED) + scalar QED with 
kinetic term of the gauge field and the Higgs degrees of freedom truncated in 
the tree level. Now that we find a gauge-invariant formulation of the Thirring 
model, we have the privilege to choose the gauge appropriate for our particular 
purpose. Here, instead of the unitary gauge notorious for loop calculations, let 
us consider the nonlocal gauge at the Lagrangian level 

where the gauge fixing parameter ^ has the momentum- (derivative-) depen- 
dence. It is straightforward to prove the possession of the BRS symmetry in 
spite of the nonlocality of ^ and thereby it guarantees the S-matrix unitar- 
ity. Another intriguing point is that, in the combined Lagrangian of Eqs. (^ 
and (||), the fictitious Nambu-Goldstone boson (j) is completely decoupled in- 
dependently of the specific form of ^(9^). Though we introduced U{1) gauge 
symmetry by use of the flavor-singlet current, the hidden local symmetry can 
easily be extended to the non-Abelian one, e.g. U (?T.)giobai x U (n)iocai symmetry 
by use of the ^V{n)/U{n)" nonlinear sigma model : 

CNon-Ab. = ^ Ipah'^D^lpa ~ ^tr [{D ■ f] , (5) 

a 

where = A'^T", and u = e'"^,!^ = ^^T", with being the U{n) gener- 
ators. Actually, Eq.(^ is gauge equivalent to the Thirring model having the 
interaction 

- ^ E i^al'-T'^^a) i^bl.T'^^b). (6) 

a, 6, a 

In contrast to the C/(l) case, however, the fictitious NG bosons in the non- 
Abelian case are not decoupled even in the gauge, which would make the 
analysis at quantum level rather complicated. 
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^Dual 



Now that we have reformulated the Thirring model as a gauge theory, we 
can further gain an insight into the theory by using a technique inherent to the 
gauge theory, namely, the dual transformation Q. We first consider the path 
integral for the Lagrangian (^), and linearize the "mass term" of gauge field 
by introducing an auxiliary field C^. Through an integration over the scalar 
field (/>,| we obtain a delta functional for S^tC^. If we pick up by use of 
dual antisymmetric-tensor field H^^...^^_^ of rank D — 2, which satisfies the 
Bianchi identity, we obtain the following path integral after integrating out the 
auxiliary field C^; 

[dH^^...^^_^][dAfj] [dV^a] [dipa] exp i J d^a;|^ ?/!aj7''^MV'o 

= j [dH,,...,,_,][di^,][dM 5C£ ^'^aT'^^^a + -^e^--''-iJ^,...^J 

expz j d^x|^Via»7^a^V^.+ ^|~^^''^^ g^,...p„_,g^" -^°-^|, (8) 

where TJ^t^ ^ ^^^ti -ff/J2 - A'-D-i ~'^M2^a'iM3 - a'-d-i "I" ' '~^(~-'^) d^^_-^H . 
The Lagrangian (||) describes N "free" fermions and a "free" antisymmetric 
tensor field of rank D ~ 2 which are, however, constrained through the delta 
functional. This implies that the dual field is actually a composite of the 
fermions. 

Let us write down the discretized Lagrange density for the Thirring model 
on the 3D Euclidean lattice. Of course, it looks not so good way to use the 
original Lagrange density with four-fermioji contact term in Eq.([^) directly. 
Therefore, one way is to discretize Eq.(^): Q 



where $ is the pseudo fermion, 6*^ is the vector auxiliary field, f3 = 1/GA is 
the inverse of the dimensionless coupling rescaled by the ultraviolet cutoff A 

■^The scalar phase <j) can in fact be divided into two parts: <j) = Q + t], where © expressed 
by multi- valued function describes the topologically nontrivial sector, e.g., the creation and 
annihilation of topological solitons, and r] given by single- valued function depicts the fluc- 
tuation around a given toppJogical sector. Inclusion of the topological sector induces a 
topological interaction term|j, though we neglect contribution in this section, since we are 
interested in <A as the Nambu-Goldstone mode. 
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which is identified with the inverse of the lattice size a, 1/A = a, and 



In Eq.(|l0|), rjfi{x) is staggered fermion phase factor, and ra a finite bare fermion 
mass term. The other method is to adopt the gauge theory formulation and 
^ke advantage of various benefits coming from the gauge theory formulation: 

Cl=Y, $^a(AftAf)-i$, - 7V/3^ cos (0(x + 0(x) + e^{x)), (11) 
where $ is the pseudo fermion and 9^ is the gauge field, 

Here 9^{x) is a gauge field which appears in the Lagrangian in terms of the 
link variable exp(i6'^). Note that the second term in Eq.(ll) is so-called the 



annealed XY model with the gauge connection between the nearest neighbors 
in statistical mechanics. 

2 Bosonization of the Thirring Model 

In (1+1) dimensions the relation in the delta functional in Eq.(||) implies 
nothing but the bosonization condition in the scheme of path integral, i.e., 
-^e'^'^d^H « -^Lip^^t^ipfj^. Integrating the fermions in Eq.(0), we obtain an 
effective theory which consists of a pseudoscalar and a vector gauge fields: 

Z2D = J [dH][dA^,] {dct ipf expi J d\ l^^d'^Hdf.H + ^^H€f,^F^"' 

(13) 

The second term of the action in Eq. (p^) is the (1+1) dimensional analogue of 
axion term which is the interaction term between the scalar and the gauge fields 
and takes the form He^ypaF^^" F^'^ in (3+1) dimensions. Though the computa- 
tion of fermionic determinant with regularization generates the Abelian chiral 
anomaly, this problem is resolved by the constant shift of scalar field H in 
axion term. Since the fermionic determinant is computed in an exact form, 
i.e., 

deti^ 1 /■ 9 ,,,, dadi,^ 

-^^i-d^=2;^y^^^^(^--i^)^' 
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the integration over gives a free massless scalar theory as the bosonized 
Thirring model 



If G is in the region — vr < G < 0, the energy per unit volume is unbounded 
below and hence the (l+l)-dimensional Thirring model with coupling constant 
G (G > or G < -tt) is well-defined. 

Though there have been several attempts to interpret the relation in Eq. 
as that for the bosonization in arbitrary dimensions a, their achievement is 
unclear yet since the fermionic determinant in spacetime dimensions more than 
(1+1)D produces nonlocal terms of the dual field. 

3 Dynamical Symmetry Breaking 

In (2+1) dimensions, our Lagrangian in Eq.(^ is invariant under the parity 



i^aix) ^ ^'^{x) = i7'7Va(x), A^{x) ^ A'^{x' ) = {- 1)' A ^ (x) . (15) 



where I]° = /, = —17^, = 7^, = —j^j^ and T" denote the genera- 
tors of U{N). The question we shall address from now on is "which symmetry is 
broken dynamically?" In concern with the parity, first issue is whether one can 
take the regularization to keep both the U{1) gauge symmetry and the parity 
or not. Since our gauge action in tree level has the parity-conserving mass and 
the number of two component Dirac fermion species is even, the parity need not 
be violated by appropriate regulator. For example, the introduction of parity- 
conserving Pauli-Villars regulator leads to the parity-invariant effective action 
for the gauge field as have done in (2+l)D quantum electrodynamics(QED3) 
Another question is whether the pattern of dynamically-generated fermion 
mass involves the parity violating mass {—Tn'ipaj^j^ipa) or not. Though at this 
stage we do not yet know whether the dynamical symmetry breaking really 
occurs or not, in this gauge- invariant formulation of the Thirring model, such 
symmetry breaking pattern is proven to be energetically unfavorable by using 
the exact argument in Ref. [10]. Namely, since the tree-level gauge action 
corresponding to Eq.(^ is real and positive semi-definite in Euclidean space, 
energetically favorable is a parity conserving configuration consisting of half 




(14) 



and the so-called global "chiral" transformation 




(16) 



5 



the 2-coinponent fcrmions acquiring equal positive masses and the other half 
equal negative masses. 

According to the above arguments, the pattern of symmetry breaking we 
shall consider is not the parity but the chiral symmetry, i.e., of which the 
breaking is U{2N) — > U{N) x U{N). Thus we investigate the dynamical mass 
of the type m4>ip in the Schwinger-Dyson equation, giving 

1 f d^q A{^q^)i + B{-q^) 

r7M-727 27^ T^r 9) ^D^"'{p - q), (17) 



NJ i{2n)D '^A'{-q')q'-B\-q' 

where the full fermion propagator is written as S(p) = i[A{—p'^)f—B{—p'^)Y~^, 
and T^{p, q) and D^^{p — q) denote the full vertex function and the full gauge 
boson propagator, respectively. Task is, by employing some appropriate ap- 
proximations, to reduce Eq.(|^ to the tractable integral equation for the mass 
function M(— p^) = B{—p'^)/A{—p'^). First, we here adopt the 1/A^ expan- 
sion for T^{p, q) and D^i,{p — q) under a nonlocal gauge, i.e., they are the 
bare vertex and the one-loop vacuum polarization of massless fermion loop at 
the 1/A^ leading order. Then the Schwinger-Dyson equation ( |l7| ) becomes the 
coupled integral equations for A(— p^) and B{—p'^). They support a trivial 
solution A{—p^) = 1 and B{—p^) = at the 1/A^ leading order, however, as 
was realized in QED3, we expect to find a nonperturbative nontrivial solution 
by examining them for finite N . 

A way is, by use of the freedom of gauge choice, to require A{—p^) = 1 in a 
Schwinger-Dyson equation for A{—p'^). Then this gauge fulfills the consistency 
between the bare vertex approximation and the Ward-Takahashi identity for 
the hidden local U{1) symmetry (or the current conservation), i.e., A{0) = 1. 
The specific form of the gauge is determined by a Schwinger-Dyson integral 
equation for A{—p^), and it reduces the coupled Schwinger-Dyson equations 
into a single equation for B(—p^) which turns out to be a mass function, i.e., 
M(-p2) B(-p2): 



1 

^ ^ NJo 



q'Bjq') 
'q^ + B\q^ 



d{q^-')K{p,q;G) ,\ , (18) 



where A is ultraviolet cutoff and the kernel K{p, q; G) is given by 

K{P. q; G) = ^^_^)^^^,J^^ru.T{n^) [ dO sin^-' d{e)[D - r;(fc2)], (19) 
with k"^ = p2 _|_ g2 _ cos 9. The gauge fixing parameter ^ is a function of k'^ 
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such as 



Vik') = ^ (20) 



(D-2) 



(?1) 



where C^^ = T (2 - -j) S -j) and 2-^1(0, 5, c; z) the hypergeomet- 

ric function. Note that the kernel K{p, q; G) is positive-definite for positive 
arguments p, q, G and is symmetric under the exchange of p and q. 

Now let us show the existence of the nontrivial solution for the Schwinger- 



Dyson equation (18) when 2 < D < A. Since we have in mind the continuous 
phase transition, the solution of our interest is the nontrivial solution which 
starts to exist without gap in the vicinity of the phase transition point. Such 
a bifurcation point is identified by the existence of an infinitesimal solution 
5B{p^) around the trivial solution B{p^) ~ 0. In terms of dimensionless vari- 
ables {p = Kx^/(^-^\ 5B{p^) = KT.{x), g = G/K^-°), Eq.(^ is reduced to 
a linearized integral equation 

S(x) = 1^' dyK{x^nD-2)ynD-2).g)^^y)^ (22) 

where am — (to/A)^^^ (0 < dm < 1) is the rescaled infrared cutoff and in 
fact TO is nothing but the dynamically generated mass by the normalization 
TO — SBiw?). We can rigorously prove that, if N is equal to the maximal eigen- 
value of the kernel (^9|), then there exists a nontrivial solution T,{x). Hence, 
for a given am, each line N(g,am) on {N,g) plane depicts a line of equal 
dynamically-generated mass to = Aam^^ Therefore, the critical line is 
defined by Ncrig) — lima-^->o N{g, am) which separates the broken phase from 
the symmetric phase. It is difficult to obtain the explicit form of the critical 
line Nci{g) for arbitrary g, however we can get it in the limit of infinite four- 



fermion coupling constant, g — > 00. In this limit, the bifurcation equation (22) 
in (2+l)D is rewritten into a differential equation 

plus the infrared boundary condition I]'((7„j) = and the ultraviolet one 
[xY.'{x) + Y.{x)]x=i = 0. When iV > A^cr = 128/37r2, there is no nontriv- 



ial solution of Eq.(23) satisfying the boundary conditions, while for N < Nc- 
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the following bifurcation solutions exist: 



E(x) 



UJ 

sm<- 



In^ 



(24) 



where uj — ^ Na/N — 1, 6 = 2uj ^ arctanw and am is given by the ultraviolet 
boundary condition: 



2 V am. 
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nn, n = 1, 2, • • • 



(25) 



The solution with n = 1 is the nodeless (ground state) solution whose scaling 
behavior is read from Eg. (|25|) : 



2S 

— = e exp 
A 



27r 



^N„/N - 1 



(26) 



128/37r^ is the same as the one in 



The critical four-fermion number Nc, 
QED3 with the nonlocal gauge. 

It is turn to comment briefly on the dynamically generated mass of the 
gauge boson and the dual transformation. The vector (gauge) boson is merely 
an auxiliary field at the tree level, however it turns out to be propagating 
by obtaining the kinetic term through fermion loop effect when the fermion 
acquires the dynamical mass. In (2+1) dimensions the pole mass My of the 
dynamical gauge boson is given by 



1 

2^ 



4™^ 



2Mv 



V , -1 Mv 
tan — — 
2m 



G- 



(27) 



which always satisfies a condition that My < 2m. Furthermore, the dual gauge 
field in Eq.(0) shares exactly the same pole structure with the gauge field 
Af^ irrespectively of the phase. 
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